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Abstract 

r~| . The Kolmogorov approach to turbulence is applied to the Burgers turbulence in 

the stochastic adhesion model of large-scale structure formation. As the perturbative 
approach to this model is unreliable, here is proposed a new, non-perturbative ap- 
proach, based on a suitable formulation of Kolmogorov's scaling laws. This approach 
suggests that the power-law exponent of the matter density two-point correlation 
function is in the range 1-1.33, but it also suggests that the adhesion model neglects 
important aspects of the gravitational dynamics. 

The large-scale structure of the universe is produced by the gravitational clustering 
of an initially homogeneous matter distribution. This process can be described by the 
Newtonian equations of motion of the matter fluid, written in comoving coordinates and 
in terms of the peculiar velocity and gravitational fields [1] . These equations are nonlinear 
and, although they can be linearized to describe the early growth of small perturbations 
of the initially homogeneous distribution, the actual structure formation takes place when 
the nonlinearity plays a major role, at the aptly called nonlinear stage of gravitational 
clustering. The nonlinearity of fluid mechanics also plays a major role in the phenomenon 
of turbulence and this is the cause of the difficulty in treating this phenomenon (often 
referred to as the "unsolved problem of classical physics" ) . Since turbulence is a key actor 
in many astrophysical scenarios, it is tempting to apply methods and ideas of turbulence 
to the study of large-scale structure formation. 

An early attempt to apply Kolmogorov's scaling laws to the origin of galaxies was made 
by Weizsacker [2], but his ideas have been since mostly restricted to intragalactic turbu- 
lence and have played no role in the study of the formation of clusters and superclusters 
of galaxies or the large-scale distribution of the dark matter. However, a popular model of 
large-scale structure formation, namely, the adhesion model [3, 4], is essentially a model of 
pressure-less turbulence, namely, the type of turbulence that occurs in strongly compress- 
ible flows and is usually called Burgers turbulence. As shown in this letter, the Kolmogorov 
approach to turbulence can be applied to cosmic Burgers turbulence, employing a suitable 
formulation of the adhesion model. 
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In cosmology, scaling laws for the velocity field, such as Kolmogorov's laws, are espe- 
cially important if they can be related to scaling laws for the matter density field, because 
the positions of astronomical objects are more easily measurable than their velocities and 
indeed have been shown to follow scaling laws. The best known scaling laws in cosmology 
have been found in the distribution of galaxies [5, 6], but evidences of similar scaling laws in 
the dark-matter distribution are found in iV-body cosmological simulations. To be precise, 
these simulations show that both the dark matter and the baryonic matter form a common 
multifractal "cosmic web" structure [7]. The cosmic web is precisely the type of structure 
predicted by the adhesion model [3, 4]. 

The cosmic web consists of sheets (Zeldovich "pancakes"), filaments and nodes, so it is 
indeed multifractal, in the sense that it is formed by objects of several dimensions, namely, 
two, one and zero dimensions. However, such a distribution needs not be scale invariant. 
Nevertheless, the actual cosmic web structure is, arguably, a self-similar multifractal [7]. 
The reason why the adhesion model does not predict a distribution with definite scaling 
properties is that the Burgers equation is integrable, in sharp contrast with the Navier- 
Stokes equation of incompressible turbulence. In other words, the chaotic properties of the 
latter are not present in the former, which nicely evolves the initial conditions. Therefore, 
the natural way of producing a self-similar cosmic web structure is to use self-similar 
initial conditions, namely, an initial Gaussian velocity distribution with a power-law power 
spectrum [8]. This type of distribution evolves to a self-similar cosmic web. 

The self-similar cosmic-web solution is obtained from the exact integral of the Burgers 
equation in the zero- viscosity limit [8]. In this limit, the Burgers equation is indeed scale 
invariant, in the sense that simultaneous space and time scalings Xx and \ 1 ~ h t and the 
induced scaling \ h u leave the equation invariant. Hence, the solution fulfills the dynamical 
scaling law 

u(x, t) = ^/(i-h) t4 ( a ./ i i/(i-h) > l), (1) 

that is, the solution at any t is obtained by scaling the solution at t — 1. This scaling 
is connected with a dynamical invariant, u 1 ^ /x, which can be identified with the specific 
dissipation rate e for h — 1/3 (the Kolmogorov scaling). One can further deduce that there 
is a homogeneity scale L(t) = t 1 ^ 1 ^^ 1/(1), such that the cosmic web structure at time t 
has only formed on scales smaller than L(t) whereas the initial homogeneous distribution 
stays on larger scales. The homogeneity scale L plays a similar role to that of the integral 
scale in Navier-Stokes turbulence. 

The definition of matter density in the adhesion model is not unique, but there is an 
"analytically convenient" definition [9] that gives rise to a straightforward relation between 
velocity and density scaling laws, because the density field is expressed in terms of the 
velocity field: 

p(x) = p det [Sij - diUj(x)\ , (2) 
where po is the constant initial density. This expression simplifies to 

Sp(x) = p(x) - p = ~Po diu\x) (3) 
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in the linear regime, \5p\ <C po- In the nonlinear regime, as shock waves form, u(x) 
becomes discontinuous and d(Uj — > — oo. The shock waves are also caustics [3, 4], where 
matter accumulates and p — > oo, according to Eq. (2). 

Unfortunately, the above-explained approach has obvious shortcomings: (i) the initial 
power spectrum is not a power law; (ii) the adhesion model is just a simplified model 
of the gravitational dynamics, which actually is chaotic; and (iii) the adhesion of matter 
into caustics is considered as an inelastic collision but the dissipated energy is not ac- 
counted for. In fact, the dissipation in the gravitational dynamics is linked to its chaotic 
nature: self-gravitating systems tend to virial equilibrium, which is independent of the 
initial conditions and implies that entropy grows in the process (this process is usually 
called virialization). Therefore, it is reasonable to supplement the adhesion model with a 
"noise", which reverts the lost kinetic energy, on the one hand, and makes the long-time 
evolution of the velocities independent of the initial conditions, on the other hand. This 
stochastic adhesion model possesses an attractor characterized by a dynamical scaling that 
is independent of the initial conditions, unlike the one defined by Eq. (1). The resulting 
stationary state, in which the energy injected on scales > L is dissipated at constant rate 
e at the Kolmogorov scale, is analogous to the stationary state of incompressible turbu- 
lence. However, in Burgers turbulence, the dissipation takes place in caustics and has 
more spatial variation than in incompressible turbulence, producing strong intermittency. 
Intermittency leads to deviations from Kolmogorov's scaling for higher order correlation 
functions. A beautiful exposition of Kolmogorov's ideas and of intermittency is given by 
Frisch [10]. 

The stochastic Burgers equation is well studied, since it appears, in terms of the ve- 
locity potential, as an equation for surface growth (the surface's height is given by the 
potential). The equation for the potential, called the KPZ equation, includes a Gaussian 
noise with power spectrum D(k,u). It has been studied with renormalized perturbation 
theory [11]. With this method, the types of noise that give rise to dynamical scaling are 
determined as fixed points of the dynamical renormalization group. For white noise, in 
three spatial dimensions, the nontrivial fixed point is repulsive, so the nonlinear term of the 
Burgers equation is irrelevant (in the renormalization group sense) and the viscous term 
dominates in the perturbative stationary state. Therefore, turbulence can only occur in the 
strong-coupling, non-perturbative regime. The addition of "colored" noise with power-law 
spectrum D(k) = Dk~ 2p [11] does not improve the situation: there can be several fixed 
points, but, in three dimensions, only the trivial fixed point is stable and only if p is small; 
otherwise, it becomes a saddle point. This means that a noise with sufficient power on 
large scales inevitably leads to a strong-coupling stationary state. 

The method of Medina et al [11] has been adapted to the cosmology setting by Dominguez 
et al [12]. They consider noise with power-law spatial and temporal correlations, D(k,u) = 
Dk~ 2p u~ 26 , and one more coupling (a sort of QFT "mass"). The corresponding renormal- 
ization group equations have several fixed points, but only one is stable. A choice of p and 
9 in certain ranges yields exponents for the power-law velocity-potential correlation func- 
tion such that the corresponding exponents of the density correlation function, obtained 
through Eq. (3), fit the range of measured exponents 7 of the two-point correlation function 
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of galaxies. This intriguing derivation of 7 has several questionable aspects, besides the 
ad-hoc choice of p and 9. First, Galilean invariance is broken, as the existence of the stable 
fixed point demands a non-zero 9 [12]. Second, Eq. (3) is only valid in the linear regime, in 
principle. Third, the results, apart from the values of 7, are also questionable: Regarding 
the values of the couplings at the fixed-point, the strength D of the correlated part of the 
noise (dominant for small k) turns out to be negative. Furthermore, the "mass" scale is 
non-vanishing, so the stable fixed point does not seem to correspond to a scale invariant 
stationary state. 

At any rate, one can argue, on general grounds, that perturbation theory (especially, the 
one-loop approximation) is not the right approach to Burgers turbulence. The effective 
coupling constants in the renormalization group equations have the generic expression 
\ 2 D/u 3 (except the "mass"), where A is the nonlinear coupling constant (to be set to the 
value of unity), D is a noise strength, and v the viscosity. Therefore, as the nonlinear 
term dominates in the inertial range, the coupling must be strong. More precisely, the 
given expression implies that the coupling constants are actually proportional to the cube 
of the Reynolds number, which has to be a very large number, making perturbation theory 
unreliable. 

Therefore, one must resort to non-perturbative methods. Standard non-perturbative 
methods in turbulence are the closure approaches, in which the hydrodynamical hierarchy 
of equations for statistical moments is closed at some order by assuming a relation between 
the moments of the corresponding order and lower order ones. There is a similar closure 
approach in cosmology, based on the BBGKY hierarchy [1]. This a second-order closure and 
it is consistent with a scaling ansatz for the two-point correlation functions, with just one 
power-law exponent, but this number remains undetermined, unless a connection with an 
initial power-law power spectrum of perturbations is assumed. As we avoid this connection, 
we prefer to follow the traditional non-perturbative methods in turbulence. They are based 
on reasonable assumptions, the simplest ones being Kolmogorov's universality assumptions, 
namely, homogeneity, isotropy, and scaling laws for the moments of longitudinal velocity 
increments [10]. These laws state that 



where Su = u(x + r/2) — u(x — r/2) and n G N. A general form of these scaling laws, 
suitable for introducing the effect of intermittency, is 



where g G R, and A does not depend on r. The effect of intermittency is given by the 
function ((q) as explained in the following. 

Kolmogorov's scaling laws are justified by employing the hierarchy of hydrodynamical 
equations, in particular, the second-order one, called the Karman-Howarth-Monin equation 
[10]. A version of this equation is valid for Burgers turbulence. An illuminating derivation 
of the equation has been given by Polyakov [13] in the one-dimensional case. Polyakov 
realizes that, in the equation for d t u 2 , the dissipation in the inertial range arises as a 



((5u ■ r/r) n ) oc {er) n /\ 



(4) 




(5) 
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field-theory anomaly, due to the non-differentiability of the velocity field. The form of the 
anomaly can be found by employing a point-splitting method. In the three-dimensional 
case, the calculation is more involved but it yields the simple result: 

d t u 2 (x) = -u\x)^^- + - lim ^- (Su'Su 2 ) . (6) 

The last term is the anomaly a(x), which would vanish if u(x) were differentiable. Remark- 
ably to derive Eq. (6), we do not need homogeneity or isotropy. Anyway, homogeneity and 
isotropy are part of Kolmogorov's universality assumptions and are natural in cosmology. 
From Eq. (6), one deduces that the average dissipation in the steady state is e = — (a)/2. 
This closure relation is an exact formulation of the n = 3 case of the scaling laws (4) 
for Burgers turbulence, analogous to Kolmogorov's "4/5" law of incompressible turbulence 
[10]. Therefore, in Eq. (5), £(3) = 1 (assuming that e is well defined in the limit v — > 0). 

If the probability P(5u) were Gaussian, then ((q) oc q, and necessarily ((q) = q/3, as in 
Eq. (4). However, intermittency manifests itself in a slower growth of ((q) for q > 3 [10]. In 
one dimension and with power-law correlated noise, the extent of intermittency in Burgers 
turbulence depends on the noise exponent, but it is always so strong that the maximum of 
C(q) is C = 1 [14, especially, Fig. 2]. Generalizing the results of Hayot and Jayaprakash [14] 
to three dimensions, in terms of the KPZ noise exponent p, the value p = 5/2 is such that 
the noise strength D has the dimensions of e and the (Burgers) noise correlation function 
is proportional to logr. This leads to the Kolmogorov scaling law ((q) = q/3 for q < 3 
[14]. The limit of the noise correlation function as r — > 0, equal to e, diverges for p < 5/2 
and is ill-defined for p > 5/2, in the limits v — > or L — > oo, respectively. In other words, 
e = J Q k 2 D(k) d 3 k diverges at k = oo or at k = and therefore is not universal. However, 
the r-dependent part of the noise correlation function is universal and proportional to r 2p ~ 5 
if 3/2 < p < 7/2, p 5/2. The values 5/2 < p < 7/2 correspond to large-scale forcing, 
such that ((q) = 1 if q > 3, but 2/3 < C(2) < 1. For p > 7/2, the r-dependent part of 
the noise correlation function is not universal and depends on scales > L; that is to say, 
it depends on the initial conditions. Furthermore, ((q) = 1 for q > 2, in this case. In 
this range of p, the stochastic adhesion model is presumably equivalent to the ordinary 
adhesion model with self-similar initial conditions. Therefore, the interesting range for the 
cosmic structure is p E (5/2, 7/2). Note that the noise correlation function does not need 
to be a power law: higher powers of r are not relevant in the inertial range and, besides, 
are not universal. 

Our next step is to calculate the density correlation function from Eq. (2), assuming 
Eq. (5) with ((q) = 1 for q > 3 and 2/3 < £(2) < 1. The expansion of the determinant in 
Eq. (2) yields: 

S p / Po = -diu 1 + (d 1 u 1 d 2 u 2 - d 1 u 2 d 2 u 1 + d^d^u 3 - d^d^u 1 + d 2 u 2 d 3 u 3 - d 2 u 3 d 3 u 2 ) 

+ 0(du) 3 . 

While Eq. (3), of 0(du), is valid in the linear regime, we have to consider the formation of 
caustics. Caustics are due to the blowing up of the eigenvalues of the matrix d{V? = dij<j>, 
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where is the velocity potential (it is discontinuous). If only one eigenvalue diverges, 
the collapse is one-dimensional and a sheet forms, so Eq. (3) is justified. For filaments or 
nodes, more terms are necessary. 

The reduced two-point correlation function of the density is 

(6p(r) 5p{Q))/pl = (d lU \r) d jU ^0)) - 2 c(r) + ( (^) 4 > + • • • , (7) 

where 

c(r) = ((d 1 u 1 d 2 u 2 - d l u 2 d 2 u 1 + d^dsu 3 - d^d^u 1 + d 2 u 2 d 3 u 3 - d 2 u 3 d 3 u 2 ){r) d jU j (0)}. 

We have shown explicitly only terms up to 0{du) 3 , because the other terms do not require 
any calculation, as we now explain. The functions on the right-hand side of Eq. (7) are 
power-laws of r, each one with a characteristic exponent —7 that can be deduced from 
Eq. (5); namely, —7 = ((n) — n for (0(du) n ). Given that ((n) — 1 for n — 3,4,5,6, we 
have 7 = 2,3,4,5, respectively. However, the maximal value is 7 = 3, which is the value 
for a Poisson distribution (shot-noise) term: this term can appear as either S(r) or r~ 3 
(see, e.g., [15]). As 7 = 3 is reached for n > 4, we only need to consider the cases n = 3 
and n = 2. 

To explicitly calculate c(r) oc r~ 2 , it is useful to express it as c(r) = Ag(r), where 

g (r) = ((<9n0 d 22 (j) - d 12 (f)d 21 (f) + d U (f)d 33 (f) - d n <f)d 31 (f) + d 22 <pd 33 <p - d 23 (f)d 32 (f))(r) 0(0)). 

This function is a dimensionless scalar, so it must be a constant. Therefore, c(r) = and 
the 7 = 2 contribution vanishes. 

In the end, the relevant contribution to the density two-point correlation function is due 
to the velocity two-point correlation function. The corresponding exponent is 7 = 2 — £(2). 
Since 2/3 < ((2) < 1, we obtain 1 < 7 < 4/3. The Kolmogorov scaling ((2) = 2/3 yields 
the upper bound, 7 = 4/3 ~ 1.33. The range of values of 7 obtained from galaxy surveys 
or iV-body cosmological simulations is (mostly) in the interval (1, 2) [5, 6, 7}. However, the 
classic value 7 = 1.7, which still stands [6, 7], is larger than 4/3. Nevertheless, a specific 
methodology for the analysis of galaxy catalogs [5] yields values of 7 in the interval (1, 1.3). 

To obtain values of 7 in the interval (4/3,2), we could take p G (3/2,5/2). Then, the 
Kolmogorov scale could not be set to zero, so it should be kept and, preferably, identified 
with a physical scale. In the gravitational dynamics, there is no intrinsic small scale, but 
there are small scales in the initial conditions. In iV-body cosmological simulations, the 
most suitable small scale is the scale of gravitational smoothing. In any case, if we were to 
take p G (3/2,5/2), then £(3) < 1, so c(r) would not vanish. When the density correlation 
function is the sum of different powers of r, the most singular term dominates in the 
nonlinear domain r ^ L (the inertial range). The most singular term is, of course, the 
Poisson term, but it must be discarded [15]. The next singular component is c(r), although 
it vanishes for p > 5/2. In contrast, for p < 5/2, c(r) would not vanish and would lead to 
a 7 > 2, instead of a 7 G (4/3, 2). 

With 1 < 7 < 4/3, the density correlation function is dominated by sheets. Interest- 
ingly, the analysis of iV-body cosmological simulations leads to a similar conclusion: the 
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bulk of mass belongs to sheets [7] . However, to fully understand the role of the three types 
of cosmic-web singularities, namely, sheets, filaments and nodes, one must go beyond the 
scope of the adhesion model, because the three types of singularities are very different 
in regard to the gravitational dynamics. The accumulation of matter in sheets leads to 
density singularities, but the gravitational potential stays finite. In contrast, filaments 
and nodes are gravitational singularities as well as density singularities, so their formation 
involves the dissipation of an infinite amount of energy. Therefore, it is not surprising 
that the analysis of iV-body cosmological simulations [7] shows that the spectrum of local 
dimensions a is cut off at a = 1, which is precisely the local dimension of filaments. For 
filaments, the gravitational potential has just a logarithmic singularity, which is milder 
than the r" 1 singularity of nodes and, hence, involves less dissipation. 

At any rate, gravitational singularities cannot be described in a Newtonian framework 
and need the Theory of General Relativity. In this theory, the energy dissipated in, for 
example, the formation of a point singularity as a black hole is finite, namely, it is given 
by the well-studied black-hole entropy. Plausibly, a good part of the gravitational energy 
dissipation at the (cosmic) Kolmogorov scale can be attributed to the formation and growth 
of super-massive black holes, which occur due to dissipative processes in the dark matter 
and, preferentially, in the baryonic matter. However, the formation and growth of black 
holes or other relativistic gravitational singularities is beyond the scope of the adhesion 
model and even beyond the scope of (state-of-the-art) iV-body cosmological simulations. 

In conclusion, a hydrodynamic closure approach to three-dimensional Burgers turbu- 
lence leads to Kolmogorov's scaling laws, although in a general form compatible with the 
presence of intermittency. These scaling laws can be applied to the stochastic adhesion 
model of the cosmic structure, in particular, to the determination of the density two-point 
correlation function. The result is in partial agreement with the two-point correlation func- 
tion obtained from the distribution of galaxies and from iV-body simulations but suggests 
that the adhesion model underestimates the contribution of low dimensional singularities 
(filaments and nodes) to energy dissipation, whereas iV-body simulations overestimate it. 
It is probably necessary to have a better modeling of small-scale dissipative processes, and 
this modeling may require ingredients from general relativity. 
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